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ABSTRACT 

Astronomical observations suggest that today's Universe is dominated by a dark en- 
ergy of unknown physical origin. One of the most notable consequences in many models 
is that dark energy should cause the expansion of the Universe to accelerate: but the 
expansion rate as a function of time has proven very difficult to measure directly. We 
present a new determination of the cosmic expansion history by combining distant 
supernovae observations with a geometrical analysis of large-scale galaxy clustering 
within the WiggleZ Dark Energy Survey, using the Alcock-Paczynski test to measure 
the distortion of standard spheres. Our result constitutes a robust and non-parametric 
measurement of the Hubble expansion rate as a function of time, which we measure 
with 10-15% precision in four bins within the redshift range 0.1 < z < 0.9. We demon- 
strate that the cosmic expansion is accelerating, in a manner independent of the pa- 
rameterization of the cosmological model (although assuming cosmic homogeneity in 
our data analysis). Furthermore, we find that this expansion history is consistent with 
a cosmological-constant dark energy. 
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1 INTRODUCTION 
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Observations by astronomers over the past fifteen years have 
suggested that the Universe is dominated by an unexpected 
component known as "dark energy", yet we still have no 
physical understanding of its existence or magnitude. De- 
termining the nature of dark energy is one of the most im- 
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portant challenges for contemporary cosmology because its 
presence implies that our understanding of the physics of 
the Universe is incomplete. Possible explanations are that 
the theory of General Relativity must be modified on cosmo- 
logical scales to generate an effective repulsive gravitational 
force, the Universe is filled with a diffuse material that acts 
with an effective negative pressure, or our interpretation of 
cosmological observations must be changed to reflect inho- 
mogeneities. 

One of the most important consequences of the first 
two scenarios is that the Universe underwent a transition 
from decelerating to accelerating expansion within the last 
few billion years. In this paper we use new observations 
and methods to map this expansion history in a model- 
independent, non-parametric and robust manner. We do as- 
sume cosmic homogeneity in our data analysis, and therefore 
we do not address the third scenario. 

In the standard cosmological analyses, involving 
datasets such as distant supernovae (Riess et al. 1998, Perl- 
mutter et al. 1999, Amanullah et al. 2010, Conley et al. 
2011), galaxy surveys (Giannantonio et al. 2008, Percival 
et al. 2010, Blake et al. 2011b) and the Cosmic Microwave 
Background radiation (Komatsu et al. 2011), the expansion 
rate of the Universe at different look-back times is not mea- 
sured directly. Accelerating expansion is a model- dependent 
implication of fitting cosmological data with parametric mod- 
els in which prior assumptions are made about the physical 
components of the Universe and their evolution with red- 
shift. However, the unknown physical nature of dark energy 
implies that we cannot yet be certain that any particular 
parametric model is correct. 

For example, measurements of distant supernovae pro- 
vide some of the strongest evidence for the existence of dark 
energy. The observed quantities are the apparent magni- 
tudes of the supernovae, which yield a relative luminosity 
distance as a function of redshift, i.e. Dl(z) Ho/c where 
Dl(z) is the luminosity distance, Ho is the local value of 
the Hubble expansion rate and c is the speed of light. These 
supernovae distances are well-fit by a model in which the 
physical contents of the Universe are divided into pressure- 
less matter with a current fractional energy density fi m , 
which dilutes with redshift as (1 + z) 3 , and a "cosmological 
constant" Qa whose energy density does not change with 
redshift. The best-fitting parameters assuming this model 
are O m ~ 0.27 and Qa ~ 0.73, implying a transition from 
decelerating to accelerating expansion at z « 0.6 (e.g. Con- 
ley et al. 2011). 

However, this conclusion is dependent on the assumed 
parameterized model. The acceleration of the expansion rate 
is not directly observed. When the supernovae data are sub- 
ject to model-independent non-parametric analysis, the ev- 
idence for accelerating expansion is much weaker (Wang 
& Tegmark 2005, Shapiro & Turner 2006, Shafieloo 2007, 
Sollerman et al. 2009, Shafieloo & Clarkson 2010). This is 
primarily because the acceleration rate is obtained as the 
second derivative of the noisy observed luminosity distances. 
Moreover, obtaining the cosmic expansion rate from a lu- 
minosity distance requires an additional assumption about 
spatial curvature. Although recent cosmological data is con- 
sistent with a spatially-flat Universe (Percival et al. 2010, 
Komatsu et al. 2011), this conclusion is again reached within 



the assumption of a parameterized model (unless we invoke 
a strong prior from inflationary cosmological models). 

The expansion history of the Universe as a function of 
redshift z is described by the evolution of the Hubble param- 
eter H(z) = (l + «) da/dt, where a(t) is the cosmic scale fac- 
tor at time t. In this paper we present a new determination of 
the function H(z) using a two-step approach. Firstly we ap- 
ply an Alcock-Paczynski measurement (Alcock & Paczynski 
1979) to the large-scale clustering of galaxies in the WiggleZ 
Dark Energy Survey (Drinkwater et al. 2010), thereby mea- 
suring the distortion parameter F(z) = (1 + z)Da(z)H(z)/c, 
where Da (z) is the physical angular-diameter distance. Sec- 
ondly we combine this measurement with the supernovae 
luminosity distances Dl,(z) Ho/c, using the equivalence of 
distance measurements Dl(z) = Da(z) (1 + z) 2 , to infer the 
value of H(z)/Ho- The result is an accurate non-parametric 
reconstruction of the cosmic expansion history, which is in- 
dependent of the assumed cosmological model (including 
spatial curvature). We present our results in the form of 
both individual measurements of H(z)/Ho in four redshift 
bins within the redshift range 0.1 < z < 0.9, and a non- 
parametric reconstruction of a continuous cosmic expansion 
history using a modified version of the iterative methods 
introduced by Shafieloo et al. (2006). 

Other methods for determining the Hubble expansion 
rate have been previously attempted. Firstly, H(z) may be 
inferred from the relative ages of passively-evolving galaxies 
(Stern et al. 2010). This is a promising technique although 
subject to many assumptions about the stellar populations 
and formation epochs of these galaxies. Secondly, the im- 
print of baryon acoustic oscillations in galaxy clustering has 
been applied as a standard ruler along the line-of-sight (Gaz- 
tanaga, Cabre & Hui 2009) . However, the current availability 
of large-scale galaxy survey data restricts this latter tech- 
nique to low redshift z — 0.35 and it is debatable whether 
the signal-to-noise of current measurements is in fact suffi- 
cient to ensure a robust measurement (Kazin et al. 2009). Fi- 
nally, the physical presence of dark energy has been inferred 
from the late-time Integrated Sachs- Wolfe effect which cor- 
relates low-redshift galaxy overdensities with Cosmic Mi- 
crowave Background (CMB) anisotropics (Giannantonio et 
al. 2008). However, this last technique does not allow re- 
covery of the detailed expansion history and the achievable 
statistical significance is limited. 

The Alcock-Paczynski test (Alcock & Paczynski 1979) 
is a geometric probe of the cosmological model based on 
the comparison of the observed tangential and radial di- 
mensions of objects which are assumed to be isotropic in 
the correct choice of model. It can be applied to the 2- 
point statistics of galaxy clustering if redshift-space distor- 
tions, the principal additional source of anisotropy, can be 
successfully modelled (Ballinger, Peacock & Heavens 1996, 
Matsubara & Suto 1996, Matsubara 2000, Seo & Eisenstein 
2003, Simpson & Peacock 2010). This method has been pre- 
viously utilized with data from the 2-degree Field Quasar 
Survey (Outram et al. 2004), but the WiggleZ survey offers 
a superior measurement owing to the reduced importance of 
shot noise. Chuang & Wang (2011) recently performed fits to 
the tangential and radial clustering of Luminous Red Galax- 
ies within the Sloan Digital Sky Survey (SDSS) at redshift 
z — 0.35, utilizing Alcock-Paczynski information. Indeed, 
a general analysis of the tangential/radial galaxy cluster- 
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ing pattern in the presence of baryon acoustic oscillations 
demonstrates how the information may be divided into an 
overall scale distortion, quantified by a distance parameter 
D v cc (D'i/H) 1 / 3 , and a warping, quantified by the Alcock- 
Paczynski distortion factor DaH, enabling the disentangling 
of Da and H (Padmanabhan & White 2008; Taruya, Saito 
& Nishimichi 2011; Kazin, Sanchez & Blanton 2011). Such 
approaches are just becoming possible with the current gen- 
eration of large-scale galaxy surveys, and will be very pow- 
erful when applied to future datasets such as the Baryon 
Oscillation Spectroscopic Survey (BOSS, Eisenstein et al. 
2011). Another recent application of the Alcock-Paczynski 
test was presented by Marinoni & Buzzi (2010), who pre- 
sented a study of the distribution of close galaxy pairs. 

The WiggleZ Dark Energy Survey (Drinkwater et al. 
2010) at the Australian Astronomical Observatory has re- 
cently provided a new large-scale galaxy redshift survey 
dataset, allowing low-redshift cosmological measurements 
in the SDSS to be extended to higher redshifts up to 
z = 0.9. In particular, the survey has yielded a measure- 
ment of the baryon acoustic peak in the clustering pattern at 
z = 0.6 (Blake et al. 2011b). In addition, Blake et al. (2011a) 
mapped the growth of cosmic structure in this dataset across 
the redshift range 0.1 < z < 0.9 for a fixed background cos- 
mological model. The new study presented in this paper 
constitutes a generalization of that analysis incorporating 
Alcock-Paczynski distortions arising in varying cosmologi- 
cal models. 

Our paper is structured as follows: in Section 2 we in- 
troduce the WiggleZ survey dataset and power spectrum 
measurements. In Section 3 we describe our implementation 
of the Alcock-Paczynski test and our model for marginal- 
izing over redshift-space distortions. In Section 4 we com- 
bine these Alcock-Paczynski measurements with supernovae 
data to deduce the expansion history H(z). In Section 5 we 
present a non-parametric reconstruction of the continuous 
cosmic expansion history, and we summarize our findings in 
Section 6. 



2 WIGGLEZ SURVEY POWER SPECTRA 

2.1 The WiggleZ Dark Energy Survey 

The WiggleZ Dark Energy Survey (Drinkwater et al. 2010) 
at the 3.9-m Anglo-Australian Telescope has provided the 
next step in large-scale spectroscopic galaxy redshift sur- 
veys, mapping a cosmic volume ~ 1 Gpc 3 over the redshift 
range < z < 1. By covering a total of about 800 deg 2 
of sky the WiggleZ survey has mapped about 100 times 
more effective cosmic volume in the z > 0.5 Universe than 
previous galaxy redshift surveys. Target galaxies were cho- 
sen by a joint selection in UV and optical wavebands, us- 
ing observations by the Galaxy Evolution Explorer satellite 
(GALEX; Martin et al. 2005) matched with ground-based 
optical imaging from the Sloan Digital Sky Survey (York et 
al. 2000) in the Northern Galactic Cap, and from the Red- 
sequence Cluster Survey 2 (RCS2) (Gilbank et al. 2011) in 
the Southern Galactic Cap. A series of magnitude and colour 
cuts (Drinkwater et al. 2010) was used to preferentially se- 
lect high-redshift star-forming galaxies with bright emission 
lines, that were then observed using the AAOmega multi- 
object spectrograph (Sharp et al. 2006) in 1-hr exposures. 



In this study we analyzed a galaxy sample drawn from our 
final set of observations, after cuts to maximize the conti- 
guity of each survey region. The sample includes a total of 
162,323 galaxy redshifts, which we divided into four differ- 
ent redshift slices of width Az = 0.2 between z — 0.1 and 
z — 0.9. The numbers of galaxies analyzed in each redshift 
slice were N = (22072, 42160, 63737, 34354). 

2.2 Power spectrum measurements 

Our analysis is based on measurements of the power spec- 
trum amplitude of WiggleZ galaxies as a function of the 
angle to the line-of-sight 9, which we parameterize by the 
variable fi = cos 9. By comparing the amplitude of power 
spectrum modes as a function of fi we can apply an Alcock- 
Paczynski measurement as proposed by Ballinger et al. 
(1996) and Matsubara & Suto (1996) and described below in 
Section 3. We measured galaxy power spectra in six indepen- 
dent survey regions, dividing each region into four redshift 
slices. At a given redshift we fitted models to the measure- 
ments in each different survey region, convolving with the 
respective survey window function. 

In order to measure the power spectrum and calculate 
this convolution we need to define the selection function of 
the survey which describes the mean galaxy density as a 
function of position. The angular and redshift dependences 
of the selection function for each WiggleZ region were deter- 
mined using the methods described by Blake et al. (2010). 
This process models several effects including the variation 
of the GALEX target density with dust and exposure time, 
the incompleteness in the current redshift catalogue, the 
variation of that incompleteness imprinted across each 2- 
degree field by constraints on the positioning of fibres and 
throughput variations with fibre position, and the depen- 
dence of the galaxy redshift distribution on optical magni- 
tude. The galaxy power spectrum was then measured using 
the optimal-weighting scheme proposed by Feldman, Kaiser 
& Peacock (1994). Our fiducial cosmological model which 
we use to convert redshifts and angular co-ordinates to dis- 
tances is a flat ACDM model with matter density fi m = 0.27. 
We therefore assume cosmic homogeneity in our data anal- 
ysis. 

Because the WiggleZ survey consists of a series of in- 
dependent narrow-angle cones (of width £s 10 degrees on 
the sky) a "flat-sky" approximation may be used in which 
the co-ordinate x-axis is aligned with the line-of-sight at 
the centre of each region; the wavevector components (k, fx) 
are then deduced from th e three-dimen sional Fourier vector 
k = (k x , k y ,k z ) as k — \Jk% + k\ + k\ and fi — \k x \/k. The 
observed power spectra were corrected for the small effect 
of redshift blunders (Blake et al. 2010). The covariance ma- 
trix of the power spectrum measurement was deduced using 
the methodology of Feldman et al. (1994), and the convo- 
lution matrices were determined from the window function 
as described in Section 2.2 of Blake et al. (2011a). Figure 1 
displays the 2D power spectrum amplitudes as a function of 
(k,n), stacked over the survey regions at each redshift. 

The effective redshifts of the power spectrum measure- 
ments in each slice were z — (0.22,0.41,0.6,0.78), deter- 
mined by weighting each pixel in the selection function by 
its contribution to the power spectrum error (Blake et al. 
2010): 
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Figure 1. The galaxy power spectrum amplitude as a function of amplitude and angle of Fourier wavevector (k,fi), determined by 
stacking observations in different WiggleZ survey regions in four redshift slices. The contours correspond to the best-fitting non-linear 
redshift-space distortion model. We note that because of the differing degrees of convolution in each region due to the window function, a 
"de-convolution" method was used to produce this plot. Before stacking, the data points were corrected by the ratio of the unconvolved 
and convolved two-dimensional power spectra corresponding to the best-fitting model, for the purposes of this visualization. In the absence 
of redshift-space distortions, the model contours would be horizontal lines if the fiducial cosmology was equal to the true cosmology. 



where n g (x) is the galaxy number density in each grid cell 
x and P g is the characteristic galaxy power spectrum am- 
plitude, which we evaluated at a scale k = 0.1 h Mpc -1 . We 
note that our results are not sensitive to the precise value 
of the effective redshift: the measured quantity is the scale 
distortion parameter F(z c g) = (1 + z c ff)DA(z c a)H(z e g) rel- 
ative to its value in the fiducial cosmological model, rather 
than the quantity F(z e g) itself. In other words, if the effec- 
tive redshift was in error by Az c g, then the systematic error 
in F(z cft + Az cf r) / Frizes + Az e g) can be negligible even 
when the difference between F(z e g) and F(z e g + Az e g) is 
not. 

3 ALCOCK-PACZYNSKI MEASUREMENT 
3.1 The Alcock-Paczynski effect 

An Alcock-Paczynski measurement (Alcock & Paczynski 
1979) is a method for constraining the cosmological model 



by comparing the observed tangential and radial dimen- 
sions of objects which are assumed to be isotropic, i.e. pos- 
sess equal co-moving tangential and radial sizes Lq. The 
observed tangential dimension is the angular projection 
A9 = Lo/[(l + z )Da(z)]. The observed radial dimension is 
the redshift projection Az = LqH(z)/c. The intrinsic size Lo 
does not need to be known in order to recover the observable 
Az/A9 — (1 + z)Da{z)H (z) /c, which is independent of any 
assumption about spatial curvature. An Alcock-Paczynski 
measurement does not necessarily need to be applied to 
cosmological "objects" but is equally valid for an isotropic 
process such as the 2-point statistics of galaxy clustering 
(Ballinger et al. 1996, Matsubara & Suto 1996). 

The relative tangential/radial distortion depends on the 
value of F(z) — (1 + z)Da(z)H(z)/c relative to the fiducial 
model (which we take as a flat A Cold Dark Matter cosmol- 
ogy with f2 m = 0.27, although the results do not depend sig- 
nificantly on this choice). We applied the Alcock-Paczynski 
methodology to the WiggleZ survey clustering by measur- 
ing the relative amplitude of power spectrum modes as a 
function of their angle to the line-of-sight (Ballinger et al. 
1996, Matsubara & Suto 1996, Simpson & Peacock 2010), 
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Figure 2. This Figure displays the joint likelihood of the Alcock-Paczynski scale distortion parameter F(z) relative to the fiducial value 
Ffid, and the growth rate quantified by f cr$(z), obtained from fits to the 2D galaxy power spectra of the WiggleZ Dark Energy Survey 
in four redshift slices. In order to produce this Figure we marginalized over the linear galaxy bias b 2 and the pairwise velocity dispersion 
a v . There is some degeneracy between F and fag but their characteristic dependence on the angle to the line-of-sight is sufficiently 
different that both parameters may be successfully extracted. The probability density is plotted as both greyscale and contours enclosing 
68% and 95% of the total likelihood. The solid circles indicate the parameter values in our fiducial cosmological model. 



assuming the underlying isotropy of these modes in the true 
cosmological model. 

In more detail: the true values of the angular diame- 
ter distance Da(z) and Hubble expansion parameter H(z) 
may differ from our adopted fiducial values, Da(z) and 
H(z), leading to two scaling factors f± = Da/ Da and 
/ll = H j H which relate the apparent, observed tangential 
and radial wavenumber components (k' ± ,k'^) to their true 
values (k±, fey) via k' ± = f±k± and fcy = f\\k\\. The apparent 
two-dimensional galaxy power spectrum is then related to 
the underlying fiducial galaxy power spectrum (neglecting 
redshift-space distortions for the moment) by 



p«(*V) = 




(following Ballinger et al. 1996, Simpson & Peacock 2010) 

where k! = \Jk'± + > m' = k'^/k' , b is an unknown 

galaxy bias factor, and P m (k) is the underlying real-space, 
isotropic matter power spectrum. We generated this power 
spectrum shape using the CAMB software (Lewis, Challi- 
nor & Lasenby 2000) with cosmological parameters consis- 
tent with the latest observations of the Cosmic Microwave 



Background radiation (Komatsu et al. 2011): matter den- 
sity f2 m = 0.27, cosmological constant £2a = 0.73, baryon 
fraction flb/^m = 0.166, Hubble parameter h = 0.71, pri- 
mordial scalar index of fluctuations n s — 0.96 and total fluc- 
tuation amplitude as ~ 0.8. We included non-linear growth 
of structure in this model using the "halofit" prescription 
in CAMB (Smith et al. 2003). These assumptions do not 
impose a significant model dependence in our analysis, as 
discussed in Section 3.4 below. 

If the underlying power spectrum P m (fc) contains fea- 
tures at known physical scales, which we can clearly de- 
tect as a function of angle, then we can extract both distor- 
tion factors (/ii, This would be equivalent to knowing 
the absolute length scale Lo introduced above, and may be 
achieved in principle using baryon acoustic oscillations as a 
standard ruler (Seo & Eisenstein 2003, Hu & Haiman 2003, 
Glazebrook & Blake 2005). However, the 2D WiggleZ power 
spectra do not possess a sufficiently high signal-to-noise ra- 
tio to resolve the imprint of baryon oscillations and they 
do not contribute any information to the fits presented here 
(although the baryon acoustic peak may be detected in the 
ID correlation function, as described by Blake et al. 2011b). 

If the power spectrum can be well-approximated by 
a scale-free power law over the range of scales of inter- 
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est, which is a valid approximation here, then we have 
no knowledge of L and can only measure f\\/f± = 
D a (z)H(z)/Da(z)H(z) = F/Ffid, where F is the observed 
scale distortion factor and ffid is its value in the fiducial 
model. This situation corresponds to the Alcock-Paczynski 
measurement (note that any constants in front of Equation 
2 are absorbed into the unknown bias factor b). In our de- 
fault implementation for determining the best-fitting value 
of F we varied the value of f± = Faa/F and fixed /y = 1. 
We checked that our results were not significantly changed 
if we instead varied /y = F/Fm, fixing f± = 1, or varied 
both f ± = y/F Rd /F and /„ = y/F/F Sd . 



3.2 Redshift-space distortion models 

The application of the Alcock-Paczynski measurement to 
large-scale galaxy clustering is complicated by the fact that 
there is a second physical cause of anisotropy in the observed 
"redshift space": the coherent, bulk flows of galaxies toward 
clusters and superclusters that induce systematic offsets in 
galaxy redshifts and hence produce a radial (but not tan- 
gential) power spectrum distortion (Kaiser 1987, Hamilton 
1992). In other words, the redshift-space clustering pattern 
is not isotropic in the true cosmological model. 

Under a set of general assumptions, independently of 
the cosmic expansion history and the FRW metric, the 
redshift-space power spectrum of a population of galaxies 
may be written as 



P g s (fc, M ) = P gg (k) - 2 f i 2 P g6 (k)+^P e e(k), 



(3) 



where P gg (k) ^= (\6 g (k)\ 2 ), P gB (k) = (S g (k)6*(k)} and 
Pge(k) = (\8(k)\ 2 ) are the isotropic galaxy-galaxy, galaxy- 
velocity and velocity-velocity power spectra (e.g. Samushia 
et al. 2011). These equations are given in terms of S g (k), 
the Fourier transform of the galaxy overdensity field, and 
6(k), the Fourier transform of the divergence of the pecu- 
liar velocity field in units of the co-moving Hubble velocity 
(Percival & White 2009). 

We assume a linear bias factor b relating the large-scale 
galaxy overdensity S g to the underlying matter overdensity 
5 predicted by theory: 8 g = b8. In this case P gg = b 2 Pss 
and P g o = bPse- Blake et al. (2011a) have shown that Wig- 
gleZ galaxies follow a distribution closely matched to the 
underlying matter (b w 1) and that the galaxy-mass cross- 
correlation is consistent with deterministic, linear bias over 
the range of scales relevant for this analysis. 

In linear perturbation theory in an FRW metric, the 
application of the continuity equation implies that Pss = 
—fPse = f 2 Pee, where / is the growth rate of structure, 
expressible in terms of the growth factor D(a) at cosmic 
scale factor a as / = dlnD/dlna, where the growth factor 
describes the evolution of the amplitude of a single perturba- 
tion, 5(a) = D(a) 8(1). In this case we recover the large-scale 
"Kaiser limit" (Kaiser 1987): 



P g s (fc,M) =b 2 P SS (k) (1 + f3 ^ 2 f 



(4) 



where /3 = f/b. In this approximation, Equation 2 can be 
modified to include redshift-space distortions as 
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(following Ballinger et al. 1996, Simpson & Peacock 2010). 
Although there is some degeneracy between /3 and the 
Alcock-Paczynski distortion parameter F, their angular sig- 
natures in the clustering model are sufficiently different 
that both parameters may be successfully extracted after 
marginalizing over the other. 

However, we do not use the redshift-space distortion 
model of Equation 5 in our default analysis because the 
Kaiser-limit approximation is only appropriate at the largest 
scales (lowest values of k) due to the non-linear growth of 
structure at smaller scales (Jennings et al. 2011, Okumura 
& Jing 2011). As a result, fitting Equation 5 to the galaxy 
power spectra could produce a systematic error in the recov- 
ered values of F and /3; we include the equations simply to 
illustrate how F and /3 may be disentangled in model-fitting. 

The signature of redshift-space distortions has been 
modelled in more detail by many authors using various tech- 
niques including empirical fitting formulae (Hatton & Cole 
1998), numerical N-body dark matter simulations (Jennings, 
Baugh & Pascoli 2011) and perturbation theory (Crocce & 
Scoccimarro 2006, Taruya, Nishimichi & Saito 2010). A full 
comparison of WiggleZ clustering data to these models was 
presented by Blake et al. (2011a). 

In the current study, we took as our fiducial redshift- 
space distortion model the fitting formulae provided by Jen- 
nings et al. (2011), which allow us to generate the func- 
tions Pse(k) and Pee(k) in Equation 3 for given values of 
the growth rate / and the non-linear matter power spec- 
trum Pss(k) — Pm(k) (which we obtained using "halofit" in 
CAMB). The Jennings et al. (2011) formulae have been con- 
structed to be valid for a range of physical dark energy mod- 
els, and do not pre-suppose any specific continuity equation 
linking velocity and density. Unlike the Kaiser limit they do 
not assume a Friedmann-Robertson- Walker (FRW) metric, 
as discussed by Samushia et al. (2011). We consider the ef- 
fect of different choices of the redshift-space distortion model 
in Section 3.4 below. 

As a further enhancement of the redshift-space distor- 
tion model, we also considered multiplying Equation 3 by an 
empirical damping function D(k,fi), representing convolu- 
tion with uncorrelated galaxy motions on small scales. The 
standard choice for this damping function is a Lorentzian 
D(k,n) = [1 + (fc||(j„) 2 ] _1 , where a v is the pairwise velocity 
dispersion, which may be varied as a free parameter but is 
set to zero for our default model. The effects of this choice on 
the Alcock-Paczynski measurements are considered in Sec- 
tion 3.4 below. 



3.3 Model fits 

We fitted our default Jennings et al. (2011) clustering model, 
parameterized by (F,f,b 2 ), to the WiggleZ Survey power 
spectra over the range of scales k < 0.2 h Mpc -1 . Figure 
2 displays the joint likelihoods of the scale distortion pa- 
rameter and growth rate in each of the redshift slices. We 
quantify the growth rate in this Figure by fa$(z), where 
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Table 1. This Table collects together the results of the various cosmological model fits to the galaxy clustering via the Alcock-Paczynski 
(AP) measurement and supernovae (SN) data presented in this study. Results are presented for four different redshift slices. 



Statistic 


Method 


O.f < z < 0.3 


0.3 < z < 0.5 


0.5 < 2 < 0.7 


0.7 < 2 < 0.9 


Effective redshift 2 




0.22 


0.41 


0.60 


0.78 


F Gd (z) 




0.23 


0.46 


0.71 


0.97 


F(z) = (l + z)D A (z)H(z)/c 


AP fit 


0.28 ±0.04 


0.44 ±0.07 


0.68 ±0.06 


0.97 ±0.12 


/<*(*) 


AP fit 


0.53 ±0.14 


0.40 ±0.13 


0.37 ±0.08 


0.49 ±0.12 


Cross-correlation in (F, f ag) 


AP fit 


0.83 


0.94 


0.89 


0.84 


D A (z)H Q /c 


SN fit 


0.210 ±0.001 


0.376 ± 0.005 


0.526 ±0.010 


0.655 ±0.015 


D A (z)H /c 


reconstruction 


0.209 ±0.001 


0.371 ±0.003 


0.517 ±0.006 


0.652 ±0.012 


H(z)/[H (l + z)] 


AP + SN fit 


1.11 ±0.17 


0.83 ±0.13 


0.81 ±0.08 


0.83 ±0.10 


H(z)/[H (l + z)] 


reconstruction 


0.91 ±0.02 


0.88 ±0.03 


0.88 ±0.08 


0.80 ±0.08 


Om(z) 


reconstruction 


0.28 ±0.05 


0.29 ±0.07 


0.31 ±0.13 


0.22 ±0.09 


l{z) 


reconstruction 


-0.3 ± 0.2 


-0.8 ±0.9 


-1.2±2.0 


6.4 ±4.9 



cts{z) = D{z) a&(0) quantifies the normalization of the mat- 
ter power spectrum at redshift z. This is a more model- 
independent observable than the growth rate itself owing to 
the degeneracy between ag(z) and the galaxy bias b in de- 
termining the overall galaxy clustering amplitude. We note 
that there is a strong correlation between F and / ag , with 
correlation coefficients r = (0.83, 0.94, 0.89, 0.84) in the four 
redshift slices, but both parameters may be successfully ex- 
tracted. The marginalized measurements of the growth rate 
are f ag(z) = (0.53±0.14, 0.40±0.13, 0.37±0.08, 0.49±0.12). 
These measurements are consistent with those reported by 
Blake et al. (2011a) who considered fits of redshift-space 
distortion models assuming F = Fad- The resulting mea- 
surements of F, marginalizing over / and b 2 , are plotted 
in Figure 3 as the solid black data points. Our result is 
F(z) = (l+z)D A (z) H(z)/c = (0.28±0.04, 0.44±0.07, 0.68± 
0.06, 0.97±0.12) in the four redshift slices. All of the different 
measurements reported in Sections 3, 4 and 5 are collected 
together in Table 1 for convenience. 

3.4 Systematics tests 

Figure 3 presents some tests of the sensitivity of our results 
to the different assumptions in the model described above. 
Firstly, we checked that our measurement of F was not sig- 
nificantly changed if we replaced our default model with a 
range of other implementations of redshift-space distortions 
discussed in the recent literature, including generating the 
density and velocity power spectra using Renormalized Per- 
turbation Theory (Crocce & Scoccimarro 2006) and adding 
the correction terms proposed by Taruya et al. (2010). In- 
deed, using the simplication of Equation 5, where Pss(k) is 
the non-linear power spectrum obtained from CAMB, pro- 
duced an unchanged result for F (within the statistical er- 
rors), as illustrated in Figure 3 (the cyan data points). 

We also considered including and excluding the pair- 
wise velocity dispersion as a free parameter. If we fit for F 
by varying fx and fixing /m , then marginalizing over the ad- 
ditional free parameter o~„ makes no difference to our results 
(illustrated by the blue data points in Figure 3). However, 
if we fit for F by varying /m and fixing f± then the resulting 
error in F is significantly increased by marginalizing over 
a v , due to the cross-talk between k\\ and a v in the damping 
term. 



• Default o Fitting to k<0.3 h/Mpc o Fitting to k<0.1 h/Mpc 

a Fitting varying a v o Fitting with P fle .P sfl =P afi ° No-BAO polynomial model 




0.2 0.4 0.6 0.8 

Redshift z 

Figure 3. This Figure quantifies the amplitude of systematic er- 
rors in our measurements of the scale distortion parameter F(z) 
in four redshift slices relative to the fiducial value Ffi<j, marginal- 
ized over the growth rate /, galaxy bias b 2 and pairwise velocity 
dispersion a v (where appropriate). The solid black data points 
show our default measurement using the redshift-space distortion 
model provided by Jennings et al. (2011) fitting to the wavenum- 
ber range < k < fc max = 0.2 h Mpc -1 . The remaining data 
points illustrate the effects of varying these assumptions: using 
fcmax = 0.3 h Mpc -1 [red], fc max = 0.1 h Mpc -1 [green], adding 
the pairwise velocity dispersion as a free parameter [blue] , fitting 
using the large-scale Kaiser limit formula of Equation 5 [cyan], 
and using the data itself to define the real-space power spectrum 
via a polynomial fit to an angle-averaged measurement of P(k) 
[magenta] . Each subsequent data point is slightly offset in redshift 
for clarity. 



Any potential systematic error in the redshift-space dis- 
tortion model is likely to worsen at smaller scales (higher 
values of k) where the modelling becomes less accurate. 
We carefully compared different choices of the fitting range 
< k < fc max , in order to ensure that any systematic error 
in the derived distortion parameter was not significant. The 
results of some of these tests are plotted in Figure 3. Our 
fiducial choice, fc max = 0.2 h Mpc -1 , was determined by the 
consideration that the measurements of F should not dif- 
fer systematically between different implementations of the 
redshift-space distortion model. Results for fc max = 0.1 and 
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0.3 h Mpc -1 are displayed as the red and green data points 
in Figure 3. As a further test we repeated the measurements 
of F for independent scale ranges k = — 0.1, 0.1 — 0.2 and 
0.2 — 0.3 h Mpc -1 ; the results were consistent within the 
statistical errors. 

Our analysis uses an underlying isotropic matter power 
spectrum that is consistent with the latest observations of 
the Cosmic Microwave Background radiation (Komatsu et 
al. 2011). In order to show that this does not introduce any 
sensitivity to the model into our analysis we checked that 
our results were unchanged if we instead generated a func- 
tion proportional to P m (fc) using a polynomial fit to the ob- 
served, angle-averaged galaxy power spectrum in the fiducial 
model; this comparison is shown as the magenta data points 
in Figure 3. The similarity between these measurements and 
our fiducial results provides further evidence that the baryon 
acoustic oscillations (which do not appear in a smooth poly- 
nomial model) are not contributing any information to the 
Alcock-Paczynski distortion fits. 

We conclude from these tests that the systematic error 
in F induced from modelling redshift-space distortions is 
much lower than the statistical error in the measurement. 
Our results are therefore insensitive to the model adopted 
(to first order), even though a series of model assumptions 
are necessary to produce these fits. 



4 DETERMINATION OF THE COSMIC 
EXPANSION HISTORY 

We converted our Alcock-Paczynski measurements of the 
scale distortion parameter F(z) = (1+z)Da{z)H (z)/c into a 
determination of the cosmic expansion history H(z) by using 
Type la supernovae (SNe la) data to fix the cosmic distance- 
redshift relation. We used the "Union-2" compilation by 
Amanullah et al. (2010) as our supernovae dataset, obtained 
from the website http://supernova.lbl.gov/Union. This 
compilation of 557 supernovae includes data from Hamuy 
et al. (1996), Riess et al. (1999, 2007), Astier et al. (2006), 
Jha et al. (2006), Wood-Vasey et al. (2007), Holtzman et al. 
(2008), Hicken et al. (2009) and Kessler et al. (2009). 

Given that the normalization of the supernova Hub- 
ble diagram M — 51og 10 /i is treated as an unknown pa- 
rameter, the supernovae data yields the relative luminosity 
distance Dl (z) Ho /c. We determined a model-independent 
value of this quantity in each redshift slice by fitting the 
distance-redshift relation as a 3rd-order polynomial in z 
(over the redshift range < z < 0.9) and marginaliz- 
ing over the values of the polynomial coefficients. We used 
the full covariance matrix of the supernovae measurements 
including systematic errors, and checked that our results 
were not significantly changed by assuming a 2nd-order or 
4th-order polynomial instead. Our results at the four red- 
shifts z = (0.22, 0.41, 0.6, 0.78) were D A (z) H /c = (0.210 ± 
0.001, 0.376 ± 0.005, 0.526 ± 0.010, 0.655 ± 0.015), where we 
have converted the luminosity distances to angular diameter 
distances assuming Dl(z) = Da{z)(1 + z) 2 . We note that 
the inclusion of the supernovae systematics covariance ma- 
trix (compared to uncorrelated errors excluding systematics) 
increases the errors in these measurements of Da Ho jc by a 
factor of two. The supernovae data points and best-fitting 
3rd-order polynomial model are displayed in Figure 4. This 



o 




^ 0.2 0.4 0.6 0.8 



Figure 4. This Figure displays the best-fitting 3rd-order polyno- 
mial to the Union-2 compilation of supernovae data, normalized 
as a plot of angular-diameter distance versus redshift. 
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Figure 5. This Figure displays our measurement of the evolution 
of the cosmic expansion rate using Alcock-Paczynski and super- 
novae data. The expansion rate is displayed using the value of 
a/ao = H(z)/[Ho(l + z)}\ accelerating expansion implies a de- 
crease in the value of this quantity with increasing redshift. The 
black data points are obtained by combining Alcock-Paczynski 
measurements of (1+ z) D a( z ) H (z) / c in four independent redshift 
slices with supernovae distance determinations of Dl(z)Ho/c at 
these redshifts, and are independent of curvature. The thicker, 
blue data points result from applying the distance reconstruc- 
tion method of Shafieloo et al. (2006) to both the supernovae and 
Alcock-Paczynski data, producing optimal errors at both low and 
high redshift but making the additional assumptions of zero spa- 
tial curvature and that Da(z) may be expressed in terms of an 
integral over 1/H(z). Predictions are plotted for three different 
models: a fiducial ACDM model with f2 m = 0.27 (solid line), 
an Einstein de-Sitter model with £7 m = 1 (dashed line), and a 
"coasting" model where a = constant (dotted line). 

model provides a good fit to the data, with a chi-squared 
statistic of 486.6 for 519 degrees of freedom. 

Combining the Alcock-Paczynski and supernovae mea- 
surements in the four redshift slices we find that 
H(z)/[H {1 + z)] = d/a = (1.11 ± 0.17, 0.83 ± 0.13, 0.81 ± 
0.08,0.83 ± 0.10) where a = H . These results are plot- 
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ted as the black error bars in Figure 5 and compared to 
three expansion history models: a ACDM cosmology with 
matter density fi m = 0.27, an Einstcin-de-Sitter (EdS) 
model with O m = 1 and a "coasting" model for which 
da/dt = Ho at all times. We find that our measurements 
are consistent with the ACDM model. Figure 5 constitutes 
our model-independent measurement of the cosmic expan- 
sion rate da/dt, with 10-15% precision in each of the four 
redshift slices. We note that the systematic errors in the 
supernovae measurements are insignificant in the determi- 
nation of H(z)/[Ho(l + z)] because the error propagation 
is dominated by the error in the Alcock-Paczynski measure- 
ment. 

Accelerating expansion implies a decrease in the value 
of da/dt with increasing redshift. This is directly observed 
in our data: if wc fit for a constant value of a/do we find 
that the range d/do < 1, containing accelerating-expansion 
models, contains 99.86% of the probability. Hence in this as- 
sessment, accelerating-expansion models are preferred with 
a statistical significance of more than 3-a. We stress that 
this is a non-parametric result, insensitive to the redshift- 
space distortion model adopted, based solely on combining 
Alcock-Paczynski measurements of Da(z)H(z) with super- 
novae measurements of Da(z) 



5 FULL DISTANCE- REDSHIFT 
RECONSTRUCTION 

In this Section we use our datasets to map out the full cos- 
mic expansion history via the non-parametric reconstruction 
method introduced by Shafieloo et al. (2006). We note that, 
unlike the preceding analysis, this approach assumes zero 
spatial curvature and that Da (z) can be expressed in terms 
of an integral over 1/H(z) (i.e., the FRW metric). Hence we 
expect to derive tighter constraints based on more assump- 
tions. 



in terms of the supernovae data points j, is given in Shafieloo 
& Clarkson (2010): 



lnD L (z)' +L =\nD L (z) 1 + N(z) x 



E 



[\n D L ( Zj )- In D L ( Zj y] 



exp 



2A 2 



where the redshift-dependent normalization is 



N( z y 



E^ exp 

j 3 



2A 2 



(6) 



(7) 



For the inverse-variance weighting Oj we followed Shafieloo 
& Clarkson (2010) and used the error in each supernova 
luminosity distance. 

We then modified the distance reconstruction method 
to utilize the residuals with respect to the Alcock-Paczynski 
data points. We initially generated an iterative correction to 
the quantity F(z) = (1 + z)Da(z)H(z)/c using an analogous 
method to Equation 6: 



AF(z) = M(z)J2 



F{z k ) - F(z k y 



exp 



2A 2 



(8) 



where k labels the Alcock-Paczynski measurements at red- 
shifts Zk, with errors Ok in F(z), and the normalization fac- 
tor is given by 



M{z)~ 



E a 2 



exp 



(Zfc 



2A 2 



(9) 



We then converted the iteration in F(z) into an iteration in 
D L {z) = (l + z) 2 D A (z) using 



D A (z) 



i+l 



D A (zY 



E(i+^'r^A,' 



(10) 



5.1 Method for reconstructing D A (z) 

We first describe the application of the Shafieloo et al. (2006) 
distance reconstruction method to SNe la data. We then out- 
line its generalization to enable the inclusion of the Alcock- 
Paczynski measurements. 

The methodology of Shafieloo et al. (2006) is an itera- 
tive approach that deduces a distance-redshift curve Da(z) 
from an initial guess model by smoothing the residuals be- 
tween the model and measurements of D A (z) and using 
these residuals to modify the guess. Iteration proceeds until 
the chi-squared statistic between model and data is mini- 
mized. The only parameter required by the method is the 
width of the smoothing kernel in redshift, A (our fiducial 
value is A = 0.1, other choices are considered below). For 
the initial guess model to commence the iteration, we chose 
a flat ACDM model with fi m = 0.27. Our results have no 
sensitivity to the choice of initial guess (as illustrated by the 
offset between the reconstructed and fiducial distance curve 
in the top left-hand panel of Figure 6). We defined the model 
on a grid of redshift steps spaced by dz = 0.01 from z = 
to z = 1. 

The equation generating iteration i + l from iteration i, 



where Az' = 0.01 is the redshift interval of the gridding 
of the model curves. The right-hand side of Equation 10 is 
derived by expressing Da(z) = ^2 , [c/H(z')]Az' , substitut- 
ing c/H(z') = (1 + z')D A (z')/F(z Z '), and taking the deriva- 
tive with respect to F. At each step in the iteration wc 
applied the residual distance correction for both the super- 
novae and Alcock-Paczynski data points, proceeding until 
the chi-squared statistic was minimized. 

We estimated the error in the distance-redshift recon- 
struction by repeating the iterative process using bootstrap 
re-samples of the datasets. The range of resulting distance- 
redshift curves across the re-samples defines the error in the 
reconstruction. We performed extensive tests using Monte 
Carlo realizations of synthetic datasets to verify that this 
procedure produced bias-free reconstructions with reliable 
error ranges. We note that this procedure differs from that 
of Shafieloo et al. (2006) and Shafieloo & Clarkson (2010), 
who used the set of distance curves [lnDi(z) 8 ] generated 
by a single application of the iteration process of Equation 
6 to estimate the error in the reconstruction. In order to 
obtain the error in the reconstructed curves we constructed 
separate bootstrap re-samples of the supernovae and Alcock- 
Paczynski data points. 
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Figure 6. This Figure shows our non-parametric reconstruction of the cosmic expansion history using Alcock-Paczynski and supernovae 
data. The four panels of this figure display our reconstructions of the distance-rcdshift relation Dj^(z), the expansion rate a/ Ho, the 
Om(z) statistic and the deceleration parameter q(z) using our adaptation of the iterative method of Shafieloo et al. (2006) and Shafieloo 
& Clarkson (2010). The distance-redshift relation in the upper left-hand panel is divided by a fiducial model for clarity, where the model 
corresponds to a flat ACDM cosmology with O m = 0.27. This fiducial model is shown as the solid line in all panels; Einstein de-Sitter 
and coasting models are also shown defined as in Figure 5. The shaded regions illustrate the 68% confidence range of the reconstructions 
of each quantity obtained using bootstrap resamples of the data. The dark-grey regions utilize a combination of the Alcock-Paczynski 
and supernovae data and the light-grey regions are based on the supernovae data alone. The redshift smoothing scale A = 0.1 is also 
illustrated. The reconstructions in each case are terminated when the SNe-only results become very noisy; this maximum redshift reduces 
with each subsequent derivative of the distance-redshift relation [i.e. is lowest for q(z)}. 



5.2 Method for reconstructing H(z), Om(z) and 

We used the reconstructed distance-redshift curves to deter- 
mine three other diagnostics that describe the cosmic expan- 
sion history. Firstly, assuming a fiat Universe we can differ- 
entiate the reconstructed distance-redshift curves to deter- 
mine the expansion rate a /do = H(z)/[Ho{l + z)]. Secondly, 
a useful diagnostic of the expansion history that can be de- 
rived from H(z) is the "Om" statistic (Sahni, Shafieloo & 
Starobinsky 2008): 

[H(z)/H f 



Om(z) 



1 



(1 



(11) 



In a spatially-flat ACDM model this statistic is constant at 
different redshifts and equal to today's value of the mat- 
ter density parameter f2 m . In universes with different cur- 
vature, or containing dark energy with different properties 
to a cosmological constant, Om(z) would evolve with red- 



shift. Finally, by a second differentiation of the distance- 
redshift curves we can obtain the deceleration parameter 
q(z) = — 'da/a 2 : 



q(z) = (l + z) 



dH(z)/dz 
H(z) 



(12) 



The confidence ranges for the reconstruction of each quan- 
tity are established using bootstrap re-sampling. 



5.3 Results of the reconstruction 

Figure 6 displays our non-parametric reconstructions of the 
four quantities Da(z), da/dt, Om(z) and q(z). The light- 
grey regions display the 68% confidence range of reconstruc- 
tions using the supernovae data alone, and the inner dark- 
grey region derives from the combination of supernovae and 
Alcock-Paczynski measurements. 

The addition of the Alcock-Paczynski data does not 
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strongly enhance the reconstruction of the distance-redshift 
relation Da(z), which is very well- measured by supernovae 
alone. However, the direct sensitivity of F(z) to H(z) pro- 
duces dramatic improvements in the reconstruction of the 
expansion-rate history, with a factor of 3 shrinking of the 
confidence ranges in H(z) for z > 0.2. At the four redshifts 
z = (0.22, 0.41, 0.6, 0.78) the reconstruction method gives us 
H(z)/[H (l+z)] = (0.91±0.02, 0.88±0.03, 0.88±0.08, 0.80± 
0.08). These measurements are correlated, but the level of 
correlation is low because the fiducial choice of the smooth- 
ing length A = 0.1 is somewhat smaller than the bin width 
Az = 0.2. 

Our measurements establish the constancy of Om(z) 
with redshift for the first time with reasonable precision, a 
key prediction of a flat ACDM model. At the four redshifts 
listed above the reconstruction method produces Om(z) = 
(0.28 ± 0.05, 0.29 ± 0.07, 0.31 ± 0.13, 0.22 ± 0.09). If we fit 
these values for the matter density in a flat ACDM model 
we find fi m = 0.29 ± 0.03, consistent with fits to the CMB 
radiation (Komatsu et al. 2011). 

We also overplot in Figure 5 the non-parametric mea- 
surements of da/dt obtained by this reconstruction method 
in the four redshift slices. Comparison of the errors with 
the direct combination of supernovae distances and Alcock- 
Paczynski data shows that at low redshift the supernovae 
data alone place the most powerful constraint on the ex- 
pansion history, greatly improving the results presented in 
Section 4, but at higher redshifts the Alcock-Paczynski data 
are crucial in order to obtain accurate measurements. 

Figure 7 illustrates the effect of changing the smoothing 
length on the reconstructed expansion history a/do, consid- 
ering values A = 0.05,0.10,0.15 and 0.20. A flat ACDM 
model with fl m = 0.27, shown as the solid line, continues to 
be a good description of the expansion history regardless of 
the smoothing scale. We checked that the error in fitting f2 m 
to these data was insensitive to A when covariances between 
different redshifts were taken into account (estimating these 
covariances using the bootstrap resamples). 

6 CONCLUSIONS 

We summarize our study as follows: 

• We have performed joint fits for the Alcock-Paczynski 
scale distortion parameter F{z) — (1 + z)Da(z)H(z) jc and 
the normalized growth rate fas(z), using 2D power spec- 
tra measured from the WiggleZ Dark Energy Survey in four 
redshift bins in the range 0.1 < z < 0.9. We can separate the 
contributions of Alcock-Paczynski and redshift-space distor- 
tions, and our result is insensitive to the range of models we 
examined for non-linear redshift-space effects when fitted to 
the wavenumber range k < 0.2 h Mpc -1 . 

• By combining the Alcock-Paczynski fits 
with luminosity-distance measurements using Type la su- 
pernovae, we determined the cosmic expansion rate d/do = 
H(z)/[Ho(l + z)\ in four redshift slices with 10-15% accu- 
racy. Our results for redshifts z = (0.22,0.41,0.6,0.78) arc 
d/d = (1-11 ± 0.17,0.83 ± 0.13,0.81 ± 0.08,0.83 ± 0.10). 
These measurements are independent of spatial curvature. 

• Our measurements show that the value of d was lower at 
high redshift, demonstrating that the cosmic expansion has 
accelerated. If we fit for a constant value of d in the redshift 



range 0.1 < z < 0.9 we find that accelerating-expansion 
models are favoured with a statistical significance of more 
than 3-cr. 

• We used an adapted version of the reconstruction 
method of Shafieloo et al. (2006) to model the continu- 
ous cosmic expansion history since z = 0.9 using both 
the Alcock-Paczynski and supernovae datasets. The Alcock- 
Paczynski measurements enable a much more accurate de- 
termination of the expansion history than the supernovae 
data alone. We demonstrate that the quantity Om(z) = 
{[H(z)/H Q ] 2 - 1}/[(1 + zf - 1] is constant with redshift, 
as expected in a spatially-flat ACDM model; fitting for the 
value of this quantity allows us produce an estimate of the 
matter-density parameter f2 m = 0.29 ±0.03, consistent with 
fits to the CMB radiation (Komatsu et al. 2011). 

We conclude that accelerating cosmic expansion can be 
recovered from cosmological data in a non-parametric man- 
ner and hence is a real physical phenomenon that must be 
accounted for by theory. Our results are independent of the 
cosmological model adopted, although our data analysis as- 
sumes a homogeneous Universe. The measured expansion 
history is well-fit by a cosmological constant which grows in 
relative importance with cosmic time. The combination of 
Alcock-Paczynski and supernovae data is a novel approach 
that enables direct observation of the Hubble expansion rate 
and will be strengthened in the future by the availability of 
new galaxy redshift survey data. 
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Figure 7. In this Figure we repeat the reconstruction of a/ Ho using both supernovae and Alcock-Paczynski data for different choices 
of the smoothing length A, the only free parameter in the reconstruction algorithm. The solid band shows the 68% confidence range 
spanned by the reconstructions and the solid red line is the prediction of a flat ACDM model with f2 m = 0.27. 
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